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Abstract  Acoustic  waves  are  used  in  SONAR  systems 
for  the  detection  and  classification  of  objects  which  may 
be  swimming  in  the  water  column,  proud  on  the  seafloor, 
or  buried  below  the  sediment.  In  some  applications,  such 
as  low  frequency  detection  and  classification  of  struc¬ 
turally  complex  objects  immersed  in  a  fluid,  signal  pro¬ 
cessing  techniques  must  be  aided  by  a-priori  model  based 
knowledge  of  the  target  echo.  Because  objects  of  inter¬ 
est  in  practical  applications  have  a  complex  geometry 
and  a  detailed  internal  structure,  consisting  of  elastic  ma¬ 
terials  and  fluid  partitions,  mathematical  techniques  ca¬ 
pable  of  dealing  with  generic  3-D  geometries,  and  ca¬ 
pable  of  coupling  different  physical  domains  (i.e.  the 
solid  and  the  fluid  domains)  must  be  used.  For  these  rea¬ 
sons,  the  Finite  Element  (FE)  Method  is  the  technique  of 
choice  for  many  researchers  in  the  field.  A  FEMlab  ap¬ 
plication  for  computing  the  scattering  from  undersea  tar¬ 
gets  is  presented,  and  guidelines  for  the  discretization  of 
the  computational  domain  in  the  presence  of  structure- 
borne  guided  elastic  waves  are  given.  It  is  found  that 
the  guided  elastic  waves  impose  strict  requirements  on 
the  meshing  of  the  elastic  structure  and  of  the  surround¬ 
ing  water  volume.  The  presented  results  are  relevant  also 
to  other  areas  of  acoustics,  where  one  is  interested  in 
modeling  scattering  or  radiation  from  elastic  structures 
and  within  elastic  structures.  Such  areas  include  non  de¬ 
structive  testing,  vehicle  noise  prediction  and  control, 
acoustic  transducer  design,  MEMS,  architectural  acous¬ 
tics,  and  medical  acoustics. 

Keywords  acoustics  •  structural  acoustics  •  resonances  • 
Lamb  waves 


1  Introduction 

The  modeling  of  scattering  and  radiation  from  complex 
3-D  structures  is  a  challenging  computational  task.  This 
paper  discusses  an  example  of  scattering  from  a  sub¬ 
merged  void  cylindrical  shell  with  hemispherical  end- 
caps  (see  Fig.  1 ),  which  exhibits  an  apparently  peculiar 
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convergence  behavior.  It  is  found  that  the  FE  conver¬ 
gence  behavior  described  in  Section  2  is  determined  by 
the  need  of  discretizing  short-wavelength  guided  elas¬ 
tic  waves  (called  Lamb  waves  [1],  [2],  [3]),  which  are 
present  in  elastic  structures  of  finite  thickness  [4],  [5]. 
Capturing  the  Lamb  wave  resonances  in  the  frequency 
domain  is  crucial  for  determining  the  time-domain  echo 
scattered  by  a  target. 

Section  3  presents  discretization  guidelines  which  en¬ 
sure  FE  convergence  at  the  Lamb  wave  resonances,  and 
hence  enable  the  researcher  to  synthesize  the  time  do¬ 
main  echoes  by  use  of  the  inverse  Fourier  transform.  For 
verification,  the  FEMlab  computational  results  are  com¬ 
pared  to  results  obtained  from  a  non  commercial  struc¬ 
tural  acoustics  finite-element  tool  [7],  developed  at  the 
NATO  Undersea  Research  Centre.  The  reference  tool  [7] 
is  based  on  a  hierarchic  polynomial  shape  function  de¬ 
velopment  framework  [8].  A  verification  study,  in  which 
the  reference  tool  [7]  has  been  successfully  compared 
to  other  standard  tools  used  in  the  underwater  acoustics 
community  has  been  published  in  [4].  In  particular,  the 
reference  solutions  shown  in  this  paper  have  been  ver¬ 
ified  in  [4].  To  obtain  the  results  presented  here,  both 
tools,  FEMlab  and  the  reference  tool,  have  received  the 
same  structured  mesh  as  an  input. 

Future  work  is  identified  and  conclusions  are  given 
in  Section  4.  The  results  show  that  the  FEMlab  mul¬ 
tiphysics  tool  has  the  capability  for  modeling  complex 
structural  acoustic  scattering,  propagation  and  radiation 
phenomena. 


2  Problem  description,  model  and  convergence 
issues 

Figure  1  .a  shows  the  cylindrical  target  made  of  steel  (mo¬ 
deled  as  isotropic  elastic  solid),  which  is  immersed  in 
water.  The  water  domain  surrounding  the  target  (Fig.  1  .a 
and  l.b)  is  modeled  with  the  3-D  acoustics  application 
mode  of  FEMlab,  and  the  target  is  modeled  with  the  3-D 
solid  stress-strain  application  mode  of  FEMlab.  Details 
on  how  the  physics  modes  in  the  fluid  and  solid  subdo¬ 
mains  are  set  up  can  be  found  in  the  FEMlab  Structural 
Mechanics  Module  documentation  [6].  The  steel  shell 
is  1mm  thick,  with  hemispherical  endcaps  of  outer  ra¬ 
dius  a  —  9.08cm.  The  end-to-end  length  of  the  cylinder 
is  L  —  60cm.  The  inside  of  the  cylinder  is  considered  to 
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Fig.  1  (a)  Dimensions  of  the  empty  cylindrical  shell,  (b)  Details  of  the  FEMlab  structured  mesh  in  the  vicinity  of  the  target,  with  incident 
plane  wave  denoted  by  arrow,  (c)  Full  mesh,  with  vertical  and  horizontal  circles  (labeled  “v-circle”  and  “h-circle”  respectively)  on  which 
the  acoustic  near  field  is  sampled. 


be  a  vacuum,  and  hence  the  natural  stress-free  conditions 
for  the  free  faces  of  solid  finite  elements  hold  on  the  in¬ 
terior  surface  of  the  cylinder.  A  plane  wave,  described 
by  a  known  function  pmc  =  exp(— iky),  where  k  is  the 
acoustic  wavenumber,  is  incident  on  the  cylinder  from 
broadside  along  the  direction  indicated  by  the  arrow  in 
Fig.  l.b. 

The  discretized  field  variable  is  scattered  pressure, 
Pscat,  which  is  defined  by  the  relation 

P  =  Pscat  T“ Pine  >  (D 

where  p  is  the  total  acoustic  pressure  in  the  fluid  medium. 
The  total  field  p,  as  well  as  the  incident  field  pmc,  both 
satisfy  the  Helmholtz  equation 

(y2+k>)p  =  o,k=W  (2) 

c 

(where  k  is  the  wavenumber,  /  is  the  frequency,  and  c 
is  the  soundspeed  of  water).  From  the  linearity  of  the 
Helmholtz  equation,  it  follows  that  also  pscat  satisfies 
Eq.(2),  and  thus  it  is  only  necessary  to  solve  for  psc at, 
instead  of  solving  for  pSCat  +  Pmc-  Hence,  Eq.  2  can  be 
written  with  /?scat  instead  of  p,  and  the  known  incident 
pressure  pmc  is  applied  via  the  boundary  conditions  at 
the  interface  between  the  solid  and  the  fluid.  This  ap¬ 
proach  makes  it  possible  to  obtain  generally  more  accu¬ 
rate  results  compared  to  those  obtained  by  solving  for  the 
total  held  p.  Physically,  the  scattered  held  is  the  quantity 
more  closely  related  to  the  monostatic  echo.  To  model 
the  scattered  held  in  FEMlab,  one  simply  uses  the  setup 
provided  by  the  FEMlab  acoustics  mode,  as  described  in 
[6].  The  only  difference  compared  to  the  approach  de¬ 
scribed  in  [6]  is  that  the  huid/solid  subdomain  coupling 
conditions  must  take  p-mc  and  the  normal  derivative  of 


Pine  with  respect  to  the  coupling  surface  into  account. 
Furthermore,  the  incident  held  must  be  set  to  zero  in  the 
spherical  radiation  BC’s  of  FEMlab,  which  are  applied 
on  the  spherical  boundary  of  the  computational  domain 
(Fig.  1  .c)  to  approximate  the  radiation  condition  for  psca t. 
The  radiation  BC  boundary  is  located  approximately  at 
a  distance  of  4  acoustic  wavelengths  from  the  cylindrical 
target. 

Because  the  problem  is  symmetric  with  respect  to 
two  planes  passing  through  the  origin  and  containing  the 
wavenumber  vector  of  the  incident  plane  wave,  with  one 
such  plane  being  perpendicular  to  the  cylinder  axis,  and 
the  other  one  being  perpendicular  to  the  cylinder  radius, 
it  is  sufficient  to  model  only  one  quarter  of  the  geome¬ 
try,  as  shown  in  Fig.  l.c.  The  symmetry  of  the  solution 
is  enforced  using  homogeneous  Neumann  BC’s  on  the 
intersections  of  the  symmetry  planes  with  the  spherical 
computational  fluid  domain,  and  using  zero  normal  elas¬ 
tic  displacement  conditions  (Dirichlet  type)  on  the  two 
surface  strips  of  finite  thickness  obtained  by  intersecting 
the  cylindrical  target  (elastic  domain)  with  the  symmetry 
planes. 

Computations  are  performed  at  three  different  fre¬ 
quencies,  namely  f\  —  2.5kHz,  —  2.895kHz  and  fo  — 

3.5kHz.  The  frequency  /2  corresponds  to  the  location 
of  a  Lamb  wave  mode  resonance  [4].  Examples  of  the 
FEMlab  output  obtained  for  computations  at  a  frequency 
of  3.5kHz  are  presented  in  Fig.  2.  The  quantity  of  inter¬ 
est  for  the  present  study  is  the  target  strength  (TS): 

TS  =  201og10(i?  \p 

scat  I),  (3) 

where  R  is  the  observation  range,  i.e.  the  radius  of  one  or 
the  other  circular  path  described  in  Fig.  l.c.  TS  is  mea¬ 
sured  in  decibels  [dB] . 


Fig.  2  FEMlab  results  at  /  =  3.5kHz.  (a)  magnitude  of  acoustic 
pressure  in  the  near  field,  (b)  von  Mises  stress  on  the  surface  of  the 
elastic  target. 


To  make  an  accurate  comparison  between  the  differ¬ 
ent  results  possible,  the  magnitude  of  acoustic  scattered 
pressure  |/>SCat(0)|  is  sampled  on  the  two  circles  defined 
in  Fig.  l.c,  and  labeled  “v-circle”  and  “h-circle”  respec¬ 
tively,  with  0  denoting  the  angle  between  the  point  of 
impact  of  the  plane  wave  and  the  observation  point. 

The  results  obtained  from  the  mesh  of  Fig.  1  at  /  = 
/i  =  2.5kFIz  are  presented  in  Fig.  3.  Flere  and  in  what 
follows,  the  label  “COMSOL”  denotes  results  obtained 
with  FEMlab,  and  the  label  “Ref”  denotes  results  ob¬ 
tained  with  the  reference  tool  [7].  The  use  of  quadratic 
polynomial  shape  functions  is  represented  by  the  label 
“p=2”  in  the  plots,  and  the  use  of  cubic  shape  functions 
is  represented  by  the  label  “p=3”.  Figure  3  shows  that  the 
convergence  is  not  completely  reached  when  using  qua¬ 
dratic  polynomial  shape  functions.  At  the  critical  reso¬ 
nance  frequency  /  =  2.895kHz,  a  complete  lack  of  con¬ 
vergence  is  evident  in  Fig.  4.  For  p=2,  the  COMSOL  re¬ 
sult  and  the  reference  result  agree  closely,  but  they  fail 
completely  in  capturing  the  resonant  behavior  described 
by  the  p=3  curve.  As  the  frequency  is  increased  to  3.5 
kHz,  the  convergence  improves  again  (Fig.  5),  indicating 
that  the  FE  convergence  is  related  to  a  wave  phenomenon 
which  is  predominant  at  the  resonance  frequency.  Spu¬ 
rious  reflections  from  the  spherical  boundary  on  which 


Fig.  3  Results  on  the  v-circle  at  /  =  2.5kHz.  Agreement  between 
“COMSOL  p=2”  (solid)  result  and  the  “Ref  p=2”  (crosses)  is  so 
close  that  the  two  curves  cannot  be  distinguished  on  the  plot. 


Fig.  4  Results  on  the  v-circle  at  resonance  /  =  2.895kHz.  The 
p=2  curves  appear  to  miss  completely  the  resonant  behavior  of  the 
structure. 


the  radiation  BC’s  are  applied  are  not  evident  in  the  re¬ 
sults,  as  can  be  seen  from  the  curves  at  2.5  kHz  and  at 
3.5  kHz,  which  suggests  that  the  radiation  conditions  are 
performing  well. 

At  a  first  glance,  the  behavior  of  the  FE  solutions  ap¬ 
pears  to  be  rather  peculiar,  since  the  elements  used  in  the 
mesh  discretize  the  surface  of  the  cylinder  with  as  many 
as  18  elements  per  acoustic  wavelength  Awatel-  =  c/f, 
with  c  =  1500m/s  being  the  soundpeed  in  the  surround¬ 
ing  water.  This  appears  to  be  in  contrast  with  the  com¬ 
mon  practice  of  discretizing  the  domain  with  6  quadra¬ 
tic  (p=2)  elements  per  Awater  (see  for  example  [9]).  Such 
meshing  guidelines  are  justified  by  the  analysis  of  the  ap¬ 
proximation  of  sinusoids  of  a  given  wavelength  using  the 
polynomial  shape  functions.  As  shown  below,  the  guide¬ 
line  is  valid  also  in  the  case  considered  here,  provided 
that,  instead  of  limiting  the  attention  only  to  Awatcl-  and 
to  evanescent  waves  which  decay  rapidly  with  distance 
from  the  target,  one  identifies  the  wavelength  of  all  the 
relevant  wave  phenomena  occurring  in  the  fluid  loaded 
structure  of  finite  thickness. 


Fig.  5  Results  on  the  v-circle  at  frequency  /  =  3.5kHz.  The  dis¬ 
agreement  between  the  p=2  curves  and  the  p=3  curve  is  not  as 
evident  as  in  the  2.895  kHz  resonance  case. 


3  Lamb  wave  based  discretization 

The  convergence  behavior  of  the  FE  model  can  be  ex¬ 
plained  by  noticing  that  the  shortest  wavelength  in  the 
problem  is  not  the  wavelength  Awatcr  —  c/ f  computed 
with  respect  to  the  soundspeed  in  the  water  c  —  1500  m/s. 
An  analysis  of  the  elastic  waves  present  in  structures  of 
finite  thickness  according  to  [2]  reveals  that  the  smallest 
wavelength  encountered  in  the  problem  is  the  one  as¬ 
sociated  with  the  asymmetric  subsonic  Lamb  wave  res¬ 
onances  in  the  cylindrical  part  of  the  target.  Analytical 
methods  can  be  used  to  estimate  the  location  of  the  Lamb 
wave  resonances  and  the  corresponding  phase  speed  of 
the  resonant  mode  for  an  ideally  infinite  cylindrical  shell, 
and  for  a  spherical  shell.  Figure  6  shows  the  analytically 
estimated  Lamb  wave  resonance  frequencies  and  associ¬ 
ated  Lamb  wave  wavelengths  for  an  infinite  steel  shell 
and  for  a  spherical  steel  shell  having  the  same  thick¬ 
ness  and  the  same  radius  as  the  cylinder  considered  here. 
While  it  can  be  expected  that  the  estimates  of  the  Lamb 
wave  wavelengths  obtained  from  Fig.  6  are  applicable 
to  the  finite  cylinder  with  hemispherical  endcaps  used  in 
the  FE  computations,  the  resonance  frequency  locations 
cannot  be  expected  to  coincide  with  the  ones  computed 
for  the  ideally  infinite  cylinder  and  for  the  sphere.  Thus, 
from  Fig.  6  it  can  be  seen  that  the  wavelength  of  the 
Lamb  wave  in  the  cylindrical  body,  (Aao-  =  0.05  m  at 
3.5kFIz)  is  approximately  8.5  times  shorter  than  Awater- 
The  phase  speed  corresponding  to  the  shortest  subsonic 
Lamb  wave  is  approximately  175  m/s,  while  the  sound- 
speed  in  water  is  1500  m/s. 

It  should  be  noted,  that  for  complex  structures  one 
cannot  determine  a  priori  whether  a  given  Lamb  wave 
mode  resonance  decays  rapidly  with  distance  from  the 
target  in  the  near  field,  or  whether  it  propagates  out  into 
the  far  field.  Experimental  and  modeling  results  show¬ 
ing  the  propagation  of  Lamb  waves  in  the  far  field  are 
presented  for  example  in  [10].  The  results  of  reference 
[10]  imply  that  the  Lamb  waves  in  the  elastic  structure 
cannot  be  neglected,  even  though  one  may  be  primarily 
interested  in  computing  only  far  field  quantities. 


Frequency  [kHz] 


Fig.  6  Wavelengths  of  subsonic  asymmetric  Aq  Lamb  waves  with 
modal  resonance  orders  n  in  the  frequency  band  of  interest  (area 
shaded  in  grey)  for  a  spherical  shell  (diamonds),  and  for  a  cylin¬ 
der  (asterisks).  Dashed  line  shows  wavelength  w.r.t.  soundspeed  in 
water. 


Since  the  Lamb  waves  travel  along  the  thin  elastic 
structure,  it  is  necessary  to  discretize  the  directions  in¬ 
plane  with  the  target  surface  using  a  sufficient  number 
of  degrees  of  freedom  with  respect  to  Aao-  Awater- 
As  observed  above,  the  FE  mesh  used  here  has  approx¬ 
imately  a  density  of  1 8  elements  per  wavelength  on  the 
target  surface  with  respect  to  Awater-  This  implies  that  the 
Lamb  waves  are  discretized  with  only  2  elements  per 
wavelength.  The  meshing  guidelines  mentioned  above 
suggest  that  this  discretization,  in  conjunction  with  the 
use  of  quadratic  (p=2)  polynomial  shape  functions,  is  in¬ 
sufficient  to  represent  the  Lamb  waves  accurately.  In  the 
case  analyzed  here,  the  problem  can  be  overcome  eas¬ 
ily  by  raising  the  p-level  in  the  mesh  from  quadratic  to 
cubic. 

Results  on  the  “v-circle”  and  on  the  “h-circle”  for 
the  three  frequencies  of  interest,  obtained  with  the  cu¬ 
bic  polynomial  meshes  using  both  FEMlab  and  the  ref¬ 
erence  tool,  are  shown  in  Figs.  7-9.  In  all  three  cases,  the 
excellent  agreement  between  the  two  software  tools  can 
be  seen.  It  has  been  shown  in  previous  work  [4]  that  the 
solution  presented  in  Figs.  7-9  is  in  agreement  with  con¬ 
verged  solutions  obtained  using  other  alternative  codes 
for  verification.  The  most  remarkable  feature  of  the  pre¬ 
sented  solutions  is  the  improvement  obtained  at  the  res¬ 
onance  frequency,  as  shown  in  Fig.  4  and  in  Fig.  8. 

For  the  solution  of  the  FE  linear  system  in  FEMlab, 
the  SPOOLES  stationary  linear  direct  solver  with  sym¬ 
metric  option  was  used.  The  chosen  preordering  algo¬ 
rithm  was  “best  of  nested  dissection  and  multisection”. 
The  FE  linear  system  consisted  of  17986  complex  de¬ 
grees  of  freedom  (DOF)  for  the  quadratic  shape  func¬ 
tion  case  (p=2),  and  of  57352  complex  DOF  for  the  cu¬ 
bic  shape  function  case  (p=3).  The  RAM  required  by  the 
solver  was  approximately  500MB  for  the  first  case,  and 
1.3GB  for  the  second  case. 


Fig.  8  Results  on  the  v-circle  and  on  the  h-circle  at  the  Lamb  wave 
resonance  frequency  /  =  2.895kHz. 


Fig.  9  Results  on  the  v-circle  and  on  the  h-circle  at  /  =  3.5kHz. 


4  Conclusions  and  further  work 

The  work  presented  here  suggests  that  converged  results 
for  structural  acoustics  problems  are  obtained  using  me¬ 
shes  capable  of  resolving  the  short  wavelength  Lamb 
waves  in  the  elastic  structure.  For  complex  objects,  other 
than  cylinders  and  spheres  or  combinations  of  those,  one 


may  attempt  to  estimate  the  Lamb  wave  resonances  and 
wavelengths  via  analytical  models,  using  the  local  curva¬ 
ture,  material  properties  and  shell  thickness  as  inputs.  An 
analytical  study  of  this  kind  can  guide  the  researcher  in 
designing  a  convergent  finite  element  mesh.  The  conver¬ 
gence  behavior  of  the  FE  solution  at  the  resonance  fre¬ 
quency,  and  the  fact  that  the  resonances  of  fluid  loaded 
elastic  structures  cannot  be  estimated  easily  a  priori,  im¬ 
ply  that  the  convergence  of  a  steady  state  FE  computa¬ 
tion  must  be  verified  not  only  at  the  highest  frequency 
but  at  all  the  frequencies  in  the  frequency  band  of  inter¬ 
est. 

As  the  thickness  of  an  elastic  structure  increases  with 
respect  to  the  wavelength  and  with  respect  to  the  ra¬ 
dius  of  curvature,  the  asymmetric  subsonic  A(j  Lamb 
waves  transition  gradually  to  Rayleigh  waves  [11], [12], 
In  the  elastic  layer,  the  Rayleigh  waves  are  confined  to 
the  close  vicinity  of  the  interface  between  the  solid  and 
other  domains  or  layers.  From  the  fluid/elastic  interface, 
the  Rayleigh  waves  can  radiate  into  the  acoustic  far  field, 
with  the  decay  governed  by  a  problem  dependent  con¬ 
stant.  Also  the  Rayleigh  wavelength  can  be  smaller  than 
Awater,  and  hence  it  is  necessary  to  consider  such  waves 
carefully  when  generating  meshes  for  thick  structures. 

From  the  computational  point  of  view,  the  restrictive 
meshing  guidelines  imposed  by  the  elastic  wave  analysis 
put  a  heavy  burden  on  the  computational  resources,  espe¬ 
cially  if  one  considers  more  complex  problems  than  the 
one  presented  in  this  paper.  For  this  reason,  it  is  desirable 
to  limit  the  extent  of  the  fluid  domain  surrounding  the 
target  in  the  FE  part  of  the  computation  to  a  minimum. 
Far  field  quantities  in  the  free  field  can  be  obtained  by 
sampling  the  FE  results  on  a  surface  surrounding  the  tar¬ 
get,  and  by  use  of  the  Helmholtz- Kirchhoff  integral  the¬ 
orem  or  by  numerical  implementations  of  potential  the¬ 
ory  [5].  For  cases  where  one  needs  to  compute  far  field 
scattering  or  radiation  from  an  object  located  in  a  com¬ 
plex  environment,  it  is  also  possible  to  exploit  special¬ 
ized  propagation  models.  Such  generally  non  FE-based 
models  are  efficient  for  predicting  the  acoustic  propaga¬ 
tion  in  a  complex  environment,  but  they  lack  the  capa¬ 
bility  of  describing  the  scattering  from  objects.  For  this 
reason,  an  efficient  hybrid  approach  which  can  be  envi¬ 
sioned  consists  of  using  the  FE  tool  for  the  near  field 
computations  in  the  vicinity  of  the  target,  and  passing 
the  FE  results  sampled  in  the  near  field  to  the  propaga¬ 
tion  tool  of  choice,  which  can  describe  the  scattered  echo 
or  the  radiated  signal  in  the  far  field  within  the  complex 
environment  [5].  Such  an  approach  makes  it  possible  to 
solve  problems  of  practical  interest  by  exploiting  opti¬ 
mally  the  resources  of  the  FE  tool. 

This  paper  demonstrates  the  potential  of  the  FEM- 
lab  multiphysics  tool  for  modeling  the  complex  acoustic 
phenomena  which  arise  in  the  presence  of  fluid  loaded 
vibrating  or  scattering  structures.  The  wealth  of  differ¬ 
ent  solvers  offered  by  FEMlab  is  especially  attractive  for 
experimenting  with  different  solution  techniques,  both 
iterative  ones  and  direct  ones.  At  this  point,  it  should 
also  be  mentioned,  that  the  addition  of  some  features  to 
future  releases  of  FEMlab  would  make  the  tool  more  ap¬ 
plicable  to  acoustics  problems  of  interest  to  researchers 


and  engineers  working  in  the  field.  The  first  such  feature 
would  be  the  addition  of  efficient  local  numerical  radi¬ 
ation  conditions.  Such  radiation  conditions  can  reduce 
drastically  the  size  of  the  domain  modeled  with  FEMlab. 
As  a  consequence,  the  FE  problem  size  decreases  and  the 
range  of  frequencies  achievable  for  computations  with  a 
given  geometry  increases.  Candidate  techniques  for  the 
enhancement  of  numerical  radiation  BC’s  could  be  some 
local  form  of  the  higher  order  Bayliss-Turkel  radiation 
BC’s,  infinite  elements  (which  are  not  boundary  condi¬ 
tions,  but  rather  non-conventional  types  of  elements),  or 
efficient  perfectly  matched  layers  (PML's).  The  second 
issue  which  should  be  addressed  is  the  meshing  of  thin 
structures,  such  as  the  cylinder  studied  here.  Such  struc¬ 
tures  are  common  not  only  in  acoustics,  but  in  many 
areas  of  engineering  and  physics.  A  capability  for  gen¬ 
erating  structured  3-D  meshes  of  thin  objects  would  be 
the  most  desirable  option,  so  that  importing  of  meshes 
is  avoided,  and  FEMlab’s  interactive  refinement  and  p- 
enrichment  flexibility  remain  available  to  the  user. 
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